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1. Introduction 

Cosmic reionization corresponds to the transition from a fully neutral intergalactic medium 
(IGM) to an (almost) fully ionized IGM caused by the UV radiation from the first luminous objects. 
Reionization is a key benchmark in cosmic structure formation, indicating the formation of the first 
luminous objects. Reionization, and the preceding 'dark ages', remain the last of the major phases 
of cosmic evolution left to explore. 

The SKA international science advisory committee identified cosmic reionization and first 
hght as one of the key science projects for a future square kilometer array. The KSP had two parts: 

• Use the 21cm line of neutral hydrogen as the most direct probe into the evolution of the neu- 
tral intergalactic medium during cosmic reionization. Such HI 21cm studies are potentially 
the most important new window on cosmology since the discovery of the CMB. 

• Observe the gas, dust, star formation, and dynamics, of the first galaxies and AGN. Obser- 
vations at cm and mm wavelengths, provide an unobscured view of galaxy formation within 
1 Gyr of the Big Bang, and are an ideal complement to the study of stars, ionized gas, and 
AGN done using near-IR telescopes. 

I should start with a mea culpa: the original SKA KSP was misnamed as study of the 'dark 
ages'. The dark ages correspond to the period prior to reionization and first light, while the SKA 
KSP focuses on essentially cosmic reionization, which signals the end of the dark ages. This review 
will focus again on reionization, although I have included a summary of recent work suggesting 
that the dark ages themselves may be a gold mine of cosmological discoveries through the HI 21cm 
studies of the pre-galactic medium, potentially studied with a low frequency radio telescope on the 
Moon (section 3.5). 

Since this KSP was first proposed (Carilli et al. 2004), there has been an explosion in theoret- 
ical studies of the expected HI 21cm signal from reionization, as well as extensive reviews written 
on the first galaxies and the process of reionization (Furlanetto, Oh, Briggs 2007; Ellis 2007; Fan, 
Carilli, Keating 2007). In this review, I will concentrate on updating predicted signals, observa- 
tional challenges, and telescopes under construction. More on the theory of reionization can be 
found in Ferrara (this volume). 

In Section 4, I present the latest results on studies of the most distant galaxies (z > 6), and a 
look toward the promise of a 10% SKA demonstrator working at short cm wavelengths for such 
studies. I begin with a brief review of the current observational constraints on cosmic reionization. 
For more detail, see Fan, CarilU, Keating (2006). 

2. Observational constraints on cosmic reionization 

The last decade has seen the first observational evidence for cosmic reionization. The primary 
results come from the Gunn-Peterson effect, ie. Lya absorption by the neutral IGM, toward the 
most distant quasars (z ~ 6), and the large scale polarization of the CMB, corresponding to Thom- 
son scattering during reionization. These observations suggest that reionization was a complex 
process, with significant variance in both space and time, starting perhaps as high as z ~ 14, with 
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Figure 1: The neutral fraction, by volume, of the intergalactic medium as a function of redshift. Current 
constraints are shown, as well as representative models for the HI evolution. Also shown are potential causes 
for uncertainty in the measurements (adapted from Fan , Carilli, Keating 2006). 



the last vestiges of the the neutral IGM being etched-away by z 6 (see review by Fan, Carilli, 
Keating 2006). 

Figure 1 shows the current state of observational constraints on the cosmic neutral fraction as 
a function of redshift. The GP effect, and related statistical measures of eg. dark gaps in QSO 
spectra (Fan et al. 2006), suggest a qualitative change in the nature of the IGM at z 6, likely 
signifying the end of reionization. At the other extreme, the CMB large scale polarization suggests 
a significant ionization fraction extending to z ~ 11. A major task in the coming decade is to 
'connect the dots' from z ~ 6 out to z ~ 14 in this figure. 

These first measurements of cosmic reionization are truly a water-shed event. However, it 
has also become clear that these important first probes of cosmic reionization are fundamentally 
limited. Figure 1 also lists potential causes for uncertainty in the measurements (adapted from Fan , 
Carilli, Keating 2006). The CMB polarization remains only a 3a result, and it represents an integral 
measure of the Thomson scattering optical depth back to recombination, and hence can be fit by 
many different reionization scenarios. For the Gunn-Peterson effect, the IGM becomes optically 
thick to Lya absorption for a neutral fraction > 10^^, and hence the diagnostic capabilities of 
this probe effectively saturate at low neutral fractions. The GP conclusions are also depend on an 
assumed clumping factor for the IGM. 
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Figure 2: The predicted HI 21cm signals versus redshift for different HI masses, assuming a constant line 
width of 300 km s^ ' . Sensitivity of low frequency telescopes are also shown, for 1000 hour integration with 
LOFAR and the GMRT, and 100 hour with SKA. 



3. HI 21cin signals from cosmic reionization 



3.1 Large Scale Structure 

It is widely recognized that the most direct and incisive means of studying cosmic reionization 
is thiough the 21cm line of neutral Hydrogen (Furlanetto et al. 2006; Carilli 2006). The study of HI 
21cm emission from cosmic reionization entails the study of large scale structure (LSS), meaning 
HI masses > 10'^ Mq. 

Figure 2 shows the signal strength expected for a given HI mass versus redshift. A constant 
line width of 300 km s^ ^ was assumed for simplicity. The key point is that, even with future large 
area telescopes such as the SKA, HI measurements at these high redshifts will be restricted to LSS, 
not individual galaxies. Fortunately, during this epoch the entire IGM can be neutral, and the LSS 
in question is not simply mass clustering, but involves a combination of structure in cosmic density, 
neutral fraction, and HI excitation temperature. 

The optical depth on the 21cm line of neutral hydrogen is given by: 



IGkBV^^TsH [z) Ts dr 

where A is the Einstein coefficient and V21 = 1420.40575 MHz (eg. Santos et al. 2005). This 
equation shows immediately the rich physics involved in studying the HI 21cm line during reion- 
ization, with T depending on the evolution of cosmic over-densities, 5 (predominantly in the linear 
regime), the neutral fraction, x^i (ie. reionization), the HI excitation, or spin, temperature, Ts, and 
the velocity structure, ^, including the Hubble flow and peculiar velocities. 
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In the Raleigh- Jeans hmit, the observed brightness temperature (relative to the CMB) due to 
the HI 21cm hne at a frequency v = V2i/(1 is given by: 

Ts^ Zkz^,« 7(l + 5).^,(l-^)(l+^)^/^mK, (3.2) 
1+z Ts 

The conversion factor from brightness temperature to specific intensity, ly, is given by: ly = 
(hi(i+z)f '^^ ~ ^^(^ +z)~^rB Jy deg~2. These equations show that for Ts ~ Tcmb one expects 
no 21cm signal. When Ts» Tcmb, the brightness temperature becomes independent of spin tem- 
perature. When Ts << Tcmb, we expect a strong negative (ie. absorption) signal against the CMB. 

The interplay between the CMB temperature, the kinetic temperature, and the spin tempera- 
ture, coupled with radiative transfer, lead to a number of interesting physical regimes for the HI 
21cm signal (Ali 2005; Barkana & Loeb 2004; Carilli 2005; Furlanetto et al. 2006; Fan, Carilli, 
Keating 2006): 

• At z > 200 equilibrium between Tcmb, Tk, and T^ is maintained by Thomson scattering off 
residual free electrons and gas collisions. In this case Ts = Tcmb and there is no 21cm signal. 

• At z ~ 30 to 200, the gas cools adiabatically, with temperature falling as (1 +2)^, ie. faster 
than the (1-i-z) for the CMB. However, the mean density is still high enough to couple Ts and 
Tk, and the HI 21cm signal would be seen in absorption against the CMB (Sethi 2005; see 
section 3.5). 

• At z 20 to 30, collisions can no longer couple Tk to Ts, and Ts again approaches Tcmb- 
However, the Lya photons from the first luminous objects (Pop 111 stars or mini-quasars), 
may induce local couphng of Tk and Ts, thereby leading to some 21cm absorption regions 
(Cen 2006). At the same time, Xrays from these same objects could lead to local IGM 
warming above Tcmb (Chen & Miralda-Escude 2003). Hence one might expect a patch- work 
of regions with no signal, absorption, and perhaps emission, in the 21cm hne. 

• At z ^ 6 to 20 all the physical processes come to play. The IGM is being warmed by hard 

Xrays from the first galaxies and black holes (Loeb & Zaldarriaga 2004; Barkana & Loeb 
2004; Ciardi & Madau 2003), as well as by weak shocks associated with structure formation 
(Furlanetto, Zaldarriaga, & Hernquist 2005), such that Tk is likely larger than Tcmb globally 
(Furlanetto et al. 2004b). Likewise, these objects are reionizing the universe, leading to a 
fundamental topological change in the IGM, from the hnear evolution of large scale structure, 
to a bubble dominated era of HE regions (Furlanetto et al. 2005a). 

3.2 The HI 21cm signals 

Global signal: The left panel in Figure 3 shows predictions of the global (all sky) increase 
in the background temperature due to the HI 21cm line from the neutral IGM (Gnedin & Shaver 
2003). The predicted HI signal peaks at roughly 20 mK above the foreground at z ~ 10. At 
higher redshift, prior to IGM warming, but allowing for Lya emission from the first luminous 
objects, the HI is seen in absorption against the CMB. Since this is an all sky signal, the sensitivity 
of the experiment is independent of telescope collecting area, and the experiment can be done 
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Figure 3: Left: Global (all sky) HI signal from reionization (Gnedin & Shaver 2003). The shaded region 
shows the expected thermal noise in a carefully controlled experiment. Center: Predicted HI 21cm bright- 
ness temperature power spectrum (in log bins) at redshifts 8 and 12 (Mcquinn et al. 2006). The thin black 
line shows the signal when density fluctuations dominate. The dashed green line shows the predicted signal 
for Xi = 0.2 at z = 12, and x, = 0.6 at z = 8, in the (Furlanetto et al. 2004) semi-analytic model. The thick 
blue line shows the SKA sensitivity in lOOOhrs. The thick red dot-dash show the sensitivity of the pathfinder 
experiment LOFAR. The cutoff at low k is set by the primary beam. Right: The simulated SKA spectrum 
of a radio continuum source at z = 10 (Carilli et al. 2002). The straight line is the intrinsic power law 
(synchrotron) spectrum of the source. The noise curve shows the effect of the 21cm line in the neutral IGM, 
including noise expected for the SKA in a 100 hour integration. 

using small area telescopes at low frequency, with very well controlled, and calibrated, frequency 
response (Bowman, Rogers, Hewitt 2007). Note that the line signal is only ~ 10^^ that of the mean 
foreground continuum emission at ~ 150 MHz. 

Power spectra: The middle panel in Figure 3 shows the predicted power spectrum of spatial 
fluctuations in the sky brightness temperature due to the HI 21cm line (Mcquinn et al. 2006). For 
power spectral analyses the sensitivity is greatly enhanced relative to direct imaging due to the fact 
that the universe is isotropic, and hence one can average the measurements in annuli in the Fourier 
(uv) domain, ie. the statistics of fluctuations along an annulus in the uv-plane are equivalent. 
Moreover, unlike the CMB, HI line studies provide spatial and redshift information, and hence the 
power spectral analysis can be performed in three dimensions. The rms fluctuations at z = 10 peak 
at about 10 mK rms on scales i ~ 5000. 

A recent analysis by Lidz et al. (2008) shows that the MWA should be able to determine the 
HI 21cm power-spectrum over roughly a decade in wavenumber, ^ ~ 0.1 to \h Mpc^^ These 
data should be able to constrain both the amplitude and slope of the power spectrum. A key 
aspect of these measurements is determination of the redshift evolution of these quantities, with 
the amplitude of the rms power in the HI 21cm fluctuations peaking at redshift when the IGM is 
roughly 50% neutral. 

Absorption toward discrete radio sources: An interesting alternative to emission studies 
is the possibility of studying smaller scale structure in the neutral IGM by looking for HI 21cm 
absorption toward the first radio-loud objects (AGN, star forming galaxies, GRBs) (Carilli et al. 
2002). The rightpanel of Figure 3 shows the predicted HI 21cm absorption signal toward a high 



6 



Cosmic reionization and first light 



C.L. Carilli 




-2-10 1 3 
X (arcmin) 



-3-10 1 3 
X (arcmin) 



-3-10 1 3 
X (aremni) 



-2 -1 
log (5Ti (mK) 



Figure 4: The simulated HI 21cm brightness temperature distribution during reionization at z = 12, 9, 7 (left 
to right) (Zaldariagga et al. 2004) 



redshift radio source due to the 'cosmic web' prior to reionization, based on numerical simulations. 
For a source at z = 10, these simulations predict an average optical depth due to 21cm absorption 
of about 1%, corresponding to the 'radio Gunn-Peterson effect', and about five narrow (few km/s) 
absorption lines per MHz with optical depths of a few to 10%. These latter lines are equivalent to 
the Ly a forest seen after reionization (a recent treatment of this problem can be found in Furlanetto 
2006). Furlanetto & Loeb (2002) predict a similar HI 21cm absorption line density due to gas in 
minihalos as that expected for the 21cm forest. 

Tomography: Figure 4 shows the expected evolution of the HI 21cm signal during reioniza- 
tion based on numerical simulations (Zaldariagga et al. 2004; see also Zahn et al. 2007; Iliev et al. 
2006; Shapiro et al. 2006, etc.). In this simulation, the HII regions caused by galaxy formation are 
seen in the redshift range z ~ 8 to 10, reaching scales up to 2' (frequency widths ~ 0.3 MHz ~ 0.5 
Mpc physical size). These regions have (negative) brightness temperatures up to 20 mK relative 
to the mean HI signal. This corresponds to 5/xJy beam^^ in a 2' beam at 140 MHz. Only the full 
SKA will able to image such structures. 

Largest Cosmic Stromgren Spheres (CSS): While direct detection of the typical structure of 
HI and HII regions may be out of reach of the near-term EoR 21cm telescopes, there is a chance that 
even this first generation of telescopes will be able to detect the rare, very large scale HII regions 
associated with luminous quasars near the end of reionization. The expected signal is ~ 20mK 
XX/// on a scale ^ 10' to 15', with a line width of ~ 1 to 2 MHz (Wyithe, Loeb, Barnes 2005). 
This corresponds to 0.5 xx/// mJy beam^', for a 15' beam at z ~ 6 to 7, where x/// is the IGM 
neutral fraction. Figure 5 shows a simulation of the expected images from near-term reionization 
telescopes, such as the MWA, for these very large structures. The key aspect of these systems is that 
the structures are three dimensional, making them much easier to detect than continuum structures 
at a similar brightness. 

Wyithe et al. (2005) predict that, for late reionization, there should be several tens of (mostly 
fossil) large (> 4Mpc physical) CSS per 10° field of view in z ~ 6 to 8. An interesting addition to 
the QSO CSS studies is the possibility that massive galaxy formation at very high redshift is highly 
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Figure 5: The simulated HI 21cm brightness temperature distribution for the largest expected structures 
during reionization, those associated with luminous QSOs, or highly clustered massive galaxy formation 
(Wyithe, Loeb, Barnes 2005). Left shows channel images with 1.25 MHz width, and right shows the spec- 
trum at the position of the source. 

clustered, occuring in rare peaks in the cosmic density field. The total number of ionizing photons 
from star formation integrated over the lifetime of the system can match that radiated by a luminous 
QSO (Li et al. 2007; Wyithe et al. 2007), and hence generate CSS comparable to those predicted 
for the QSOs. The star forming systems will both increase the size of the spheres around bright 
QSOs (since QSO and massive galaxy formation are likely to be coupled), as well as increase the 
number of large spheres, in cases where the QSO has not turned-on yet. 

To add a personal bias to this review, my prediction is that the first robust detection of the 
neutral IGM using the 21cm line will not be through power-spectral studies, but come through a 
repeatable, very large, three dimensional 'hole' seen in images made in the very wide, deep field 
(A0 > lO**, Az > 2), with up-coming reionization telescopes. 

3.3 Telescopes 

Table 1 summarizes the current experiments under construction to study the HI 21cm signal 
from cosmic reionization. These experiments vary from single dipole antennas to study the all-sky 
signal, to 10,000 dipole arrays to perform the power spectral analysis, and potentially to image the 
largest structures during reionization (eg. the quasar Stromgren spheres). 

Most of the experiments have a few to 10% of the collecting area of the SKA, and there are 
many common features: 

• All rely on some form of a wide-band dipole or spiral antenna, eg. log-periodic yagis, sleeve 
dipoles, or bow-ties, with steering of the array response through electronic phasing of the 
elements. 

• The front-end electronics are relatively simple (uncooled amplifier/balun), since the system 
performance is dominated by the sky brightness temperature. 

• Most rely on a grouping of dipoles into 'tiles' or 'stations', to decrease the field-of-view, and 
to decrease the data rate into the correlator to a managable level. 

• The large number of array elements, and the need for wide-field, high dynamic range imaging 
over an octave, or more, of bandwidth, demand major computing resources, both for basic 
cross-correlation, and subsequent calibration, imaging, and analysis. For example, assuming 
an array of 1000 tiles, a 100 MHz bandwidth, and 8 bit samphng, the total data rate coming 
into the correlator is 1.6 Tbit s~^ The LOFAR array is working with IBM to apply the 27.4 
Tflop Blue Gene supercomputing technology to interferometric imaging (Falcke 2006). 
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Figure 6: Left: A few elements of the PAPER test array in Western Australia (Mastrantonio & Backer 2007 
in prep). Right: A full sky image fromPAPER-GB showing the brightest sources (Cas A, Cygnus A, Sun). 

First results from these path-finder telescopes are expected within the next few years. Expe- 
rience from these observations in dealing with the interference, ionosphere, and wide-field imag- 
ing/dynamic range problems will provide critical information for future experiments, such as the 
Square Kilometer Array, or a low frequency radio telescope on the moon. 

As an example. Figure 6 shows dipole elements and an all sky image from the Precision Array 
to Prope the Epoch of Reionization (PAPER). This array is being built in Western Australia, with 
a test installation in Green Bank, with the purpose of performing near-term studies of the HI 21cm 
signal from reionization. 

3.4 Observational challenges to HI 21cm studies of cosmic reionization 

Foregrounds: The HI 21cm signal from reionization must be detected on top of a much larger 
synchrotron signal from foreground emission. This foreground includes discrete radio galaxies, 
and large scale emission from our own Galaxy. The expected HI 21cm signal is about 10~^ of the 
foreground emission at 140 MHz. 

di Matteo et al. (2002) show that, even if point sources can be removed to the level of 1 /xJy, 
the rms fluctutions on spatial scales < 10' {I > 1000) due to residual radio point sources will be 
> 10 mK just due to Poisson noise, increasing by a factor 100 if the sources are strongly clustered 
(see also di Matteo, Ciardi, & Miniati 2004; Oh & Mack 2003). 

A key point is that the foreground emission should be smooth in frequency, predominantly 
the sum of power-law, or perhaps gently curving, non-thermal spectra. A number of complimen- 
tary approaches have been presented for foreground removal (Morales, Bowman, & Hewitt 2005). 
Gnedin & Shaver (2003) and Wang et al (2005) consider fitting smooth spectral models (power 
laws or low order polynomials in log space) to the observed visibilities or images. Morales & He- 
witt (2003) and Morales (2004) present a 3D Fourier analysis of the measured visibilities, where 
the third dimension is frequency. The different symmetries in this 3D space for the signal arising 
from the noise-like HI emission, versus the smooth (in frequency) foreground emission, can be a 
powerful means of differentiating between foreground emission and the EoR line signal. Santos et 
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Site 


Type 


V range 


Area 


Date 


Goal 








MHz 


m2 






Mark F 


Australia 


spiral 


100-200 


few 


2007 


All Sky 


EDGES'' 


Australia 


four-point 


100-200 


few 


2007 


All Sky 


GMRT' 


India 


parabola array 


150-165 


4e4 


2007 


CSS^ 


PAPERS' 


Australia 


dipole array 


110-190 


le3 


2008 


PS/CSS/Abs 


21CMA^ 


China 


dipole array 


70-200 


le4 


2007 


PS 


MWAd« 


Australia 


aperture array 


80-300 


le4 


2008 


PS/CSS/Abs 


LOFAR'' 


Netherlands 


aperture array 


115-240 


le5 


2008 


PS/CSS/Abs 


SKA' 


? 


aperture array 


100-200 


le6 


2015 


Imaging 



Table 1: HI 21cm Cosmic Reionization Experiments 



°http://www.atnf.csiro.au/news/newsletter/jun05/Cosmological_re-ionization.htm 

^http://www.haystack.mit.edu/ast/arrays/Edges/index.html 

'^http://gmrt.ncra.tifr.res.in/ 

^CSS = Cosmic Stromgren Spheres, PS = power spectrum, Abs = absorption 

''http://astro.berkeley.edu/'-^dbacker/EoR/ 

f http://cosmo.bao . ac.cn/proj ect.html 

^http://www.haystack.mit.edu/ast/arrays/mwa/ 

''http://www.lofar.org/ 

'http://www.skatelescope.org/ 



al. (2005), Bharadwaj & Ali (2005), Zaldariagga et al. (2004a) perform a similar analysis, only 
in the complementary Fourier space, meaning cross correlation of spectral channels. They show 
that the 21cm signal will effectively decorrelate for channel separations > 1 MHz, while the fore- 
grounds do not. The overall conclusion of these methods is that spectral decomposition should be 
adequate to separate synchrotron foregrounds from the HI 21cm signal from reionization at the mK 
level. 

Glesser, Nusser, & Bensor (2008) recently considered the relative contribution of various con- 
tributions to the mean sky brightness temperature at low frequency. They show that roughly 90% 
comes from the Galactic synchrotron emission, close to 10% from extragalactic radio sources, and 
up to ~ 2% from Galactic free-free emission. They also show that these foregrounds can be re- 
moved effectively through image reconstruction techniques that take into consideration the very 
different frequency structure of the cosmological line signal and and the foreground continuum 
emission. 

Ionosphere: A second potential challenge to low frequency imaging over wide fields is phase 
fluctuations caused by the ionosphere. These fluctuations are due to index of refraction fluctuations 
in the ionized plasma, and behave as oc v~^. 

Figure 7 shows an example of ionospheric phase errors on source positions using VLA data 
at 74 MHz (Cotton et al. 2004; Lane et al. 2004). Positions of five sources are shown for a 
series of snap shot images over 10 hours. A number of interesting phenomena can be seen. First, 
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Figure 7: The positions of five sources in the Virgo A field at 74 MHz observed with the VLA over 10 
hours (Cotton et al. 2004; Lane et al. 2004). The source positions vary with time due to fluctuations in the 
electronic path-length through the ionosphere. 



the sources are slowly moving in position over time, by ±50" over timescales of hours. These 
position shifts reflect the changing electronic path-length due to the fluctuating ionosphere (ie. tilts 
in the incoming wavefront due to propagation delay). Second, the individual sources move roughly 
independently. This is a demonstration of the 'isoplanatic patch' problem, ie. the excess electrical 
path-length is different in different directions. At 74 MHz, the typical coherent patch size is about 
3" to A" . Celestial calibrators further than this distance from a target source no longer give a vahd 
solution for the combined instrumental and propagation delay term required to image the target 
source. And third, at the end of the observation there occurs an ionospheric storm, or traveling 
ionospheric disturbance, which effectively precludes coherent imaging during the event. 

New wide field self-calibration techniques, involving multiple phase solutions over the field, 
or a 'rubber screen' phase model (Cotton et al. 2004), are being developed that should allow for 
self-calibration over wide fields. 

Interference: Perhaps the most difficult problem facing low frequency radio astronomy is 
terrestrial (man-made) interference (RFI). The relevant frequency range corresponds to 7 to 200 
MHz {z = 200 to 6). These are not protected frequency bands, and commercial allocations include 
everything from broadcast radio and television, to fixed and mobile communications. 

Many groups are pursuing methods for RFI mitigation and excision (see EUingson 2004). 
These include: (i) using a reference horn, or one beam of a phased array, for constant monitoring 
of known, strong, RFI signals, (ii) conversely, arranging interferometric phases to produce a null at 
the position of the RFI source, and (iii) real-time RFI excision using advanced filtering techniques 
in time and frequency, of digitized signals both pre- and post-correlation. The latter requires very 
high dynamic range (many bit sampling), and very high frequency and time resolution. 

In the end, the most effective means of reducing interference is to go to the remotest sites. 
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Figure 8: The interference environment in Australia. The FM band is between 87.5 and 108 MHz. The TV 
stations are between 180 and 200MHz, including broad band digital transmissions. Shown courtesy of A. 
Chippendale and R. Beresford (taken as part of the ATNF SKA Site Monitoring Program.) 

Figure 8 shows the interference measured at three different locations with population densities 
changing by roughly four orders of magnitude (Chippendale & Beresford 2006 ATNF RFI moni- 
toring program). The MWA, PAPER, and PAST have selected sites in remote regions of Western 
Australia, and China, because of known low RFI environments. Of course, the ultimate location 
would be the far-side of the moon (section 3.5). 

The technical challenges to HI 21cm observations of reionization are many. Use of spectral 
decomposition to remove the foregrounds requires careful control of the synthesize beam as a 
function of frequency, with the optimal solution being a telescope design where the synthesized 
beam is invariant as a function of frequency (Bowman et al. 2007). High dynamic range front 
ends are required to avoid saturation in cases of strong interference, while fine spectral sampling is 
required to avoid Gibbs ringing in the spectral response. The polarization response must be stable 
and well calibrated to remove polarized foregrounds (Haverkom et al. 2004). Calibation in the 
presence of a structured ionospheric phase screen requires new wide field cahbration techniques. 
The very high data rate expected for many element (> 10^) arrays requires new methods for data 
transmission, cross correlation, and storage. 

3.5 The lunatic fringe: A PAPER Moon 

At the lowest frequencies, < 50MHz or so, where we hope to study the Pre-galactic medium 
(PGM), phase fluctuations and the opacity of the ionosphere become problematic. There is long 
standing interest in building a low frequency radio telescope on the far side of the moon (Gor- 
golewski 1965; Burke 1985; Kuiper et al. 1990; Bums & Asbell 1991; Woan et al. 1997). The 
reasons are clear (section 2): no ionosphere, and sheilding from terrestrial radio frequency inter- 
ference (RFI). Two factors have acted to rekindle interest in low frequency radio astronomy from 
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Figure 9: Radio power received by the lunar orbiting radio astronomy explorer satellite in 1973 (RAE2; 
Alexander et al. 1975). The power is dominated by the Earth's auroral emission except during immersion, 
when the Earth is totally blocked by the moon. 



moon. First, scientifically, efforts to study the PGM using the redshifted HI 21cm line have spurred 
numerous ground-based low frequency projects (section 3). And second is the new NASA initiative 
to return Man to the moon, and beyond. 

While not a scientific rationale, it should be pointed out that a low frequency telescope may be 
the easiest astronomical facility to deploy and maintain on the moon. The antennas and electronics 
are high tolerance, with wavelengths > 1.5m and system noise characteristics dominated by the 
Galactic foreground radiation. Deployment could be automated, using either javelin deployment 
(EADS/ASTRON), rollout of thin polyimide fihns with metalic deposits (ROLSS; Lazio et al. 
2006), inflatable dipoles, or deployment by rovers. Likewise, being a phased array, low frequency 
telescopes are electronically steered, and hence have no moving parts. Lastly, there is no potential 
difficulty with lunar dust affecting the optics. 

Figure 9 shows radio power received by a low frequency radio receiver on a lunar orbiter from 
the 1970's (Alexander et al. 1975). The strong power received from Earth's auroral emission is 
completely gone during immersion. 

At the very low frequencies corresponding to the dark ages, going to the Moon becomes im- 
perative, due to the rapidly increasing ionospheric opacity and phase effects. A number of studies 
have considered this very low frequency, pre-reionization HI 21cm signal (Loeb & Zaldarriaga 
2004; Cen 2006; Barkana & Loeb 2005b; Shethi 2005). The HI 2Icm measurements can explore 
this physical regime at z 50 to 300, or redshifts z > 30. 

In this redshift regime the HI generally follows linear density fluctuations, and hence the ex- 
periments are as clean as CMB studies, and Tg < Tcmb, so a relatively strong absorption signal 
might be expected. Also, Silk damping, or photon diffusion, erases structures on scales i > 2000 in 
the CMB at recombination, corresponding to comoving scales = 22 Mpc, leaving the 21cm studies 
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as the only current method capable of probing to very large i in the linear regime. The predicted 
rms brightness temperature fluctuations are 1 to 10 mK on scales I = 10^ to 10^ (0.2" to 1"). These 
observations could provide the best test of non-Gaussianity of density fluctuations, and constrain 
the running power law index of mass fluctutions to large providing important tests of inflationary 
structure formation. Sethi (2005) also suggests that a large global signal, up to -0.05 K, might be 
expected for this redshift range. 

The difficulty in this case is one of sensitivity. The sky temperature is > lO^K, and using the 
equations in section 3.1, it is easy to show that, even if structures as bright at lOmK on scales of 
a few arcminutes exist, it would require '--'10 square kilometers of collecting area to detect them. 
For the more typical small scale structure being considered (ie. '--' 10"), the required collecting 
area increases to 3.6 x 10^° m^. The dynamic range requirements also become extreme, > 10^. 
It should also be kept in mind that the sky becomes highly scattered due to propagation through 
interstellar and interplanetary plasma, with source sizes obeying: Qmm ~ 1( imhz )~^'^^§- Hence, all 
objects in the sky are smeared to > 4" for frequencies < 30 MHz. 

3.6 The HI n=2 fine structure line: An another probe of the neurtal IGM? 

Recently, Sethi et al. (2007) have proposed alternative study of the HI 2^1/2 — 2^3/2 transition 
{in absorption), with 9911 MHz rest frequency.^ This transition has been notoriously difficult 
to detect in the nearby universe for two principal reasons related to the (large) Lya Einstein A 
coefficient, ^ \(f s^^ (Dennison et al. 2005). First, only ~ 10~^^ of atoms will be in the excited 
2p state, since the lifetime is only a few nanoseconds; emission of a radio photon and decay to the 
Is state is improbable. Second, the uncertainty principle dictates a ~ 100 MHz line width, due to 
the characteristic time scale for decay. Detecting such weak, broad lines is difficult. Additional 
factors relevant to the nearby universe include local dust absorption of Lya photons that would 
otherwise pump the 2p state and colhsional couphng of 2s and 2p in dense regions. 

Sethi et al. argue that the situation for the pristine neutral IGM at high redshift may be dif- 
ferent. Gunn-Peterson absorption toward quasars at z > 6 suggests that they are surrounded by 
"Cosmological Stromgren Spheres" (CSS), on physical scales of '--' 5Mpc or '--' 15' (Fan et al. 
2006). The IGM outside the CSS could still be substantially neutral and relatively unperturbed. 
Non-ionizing UV photons between Lya and the Lyman hmit would propagate into the neutral 
medium, be absorbed, and be re-emitted as hne radiation (e.g., resonant scattering of Lya; degra- 
dation of Lyj3 into Ha; two-photon decay 2s — \s). Sethi et al. show that through Une scattering 
processes a tiny population in the metastable 2s state (A '--' 8 s^^) can be maintained '^10^ times 
longer than without the photon field excited by the quasar. Unhke the 21 -cm transition, the fine 
structure transition would be seen in absorption against the microwave background. 

Sethi et al. perform a detailed calculation of the expected brightness of the 2^1/2 — '^P^ji 
transition. For an assumed UV radiation field and physical conditions in the IGM just outside a 
CSS around a z '^ 6 quasar, they predict a brightness temperature of Tg ^ —20//// /xK. The prospect 
of 991 1 MHz emission is promising for late reionization because it would be redshifted into the 1.4 
GHz band, where existing dish antennas, and optimized feed and receiver systems, are available. In 

^Transitions among hyperflne components are at 9852, 9876, and 10030 MHz in a ratio of 1:5:2; Ershov 1987. 
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contrast, study of the redshifted HI 21cm Une requires new facilities, VHF operation, and low-gain 
dipole feeds. Most importantly, the foregrounds are dramatically reduced at higher frequency. 

There are some significant uncertainties in these calculations. The Lya rate of 10^^ s^^ as- 
sumed by Sethi et al. is lOx higher than has been commonly estimated for high redshift quasars 
(e.g., Wyithe, Loeb, & Carilli 2005; Fan et al. 2006). However, re-estimation of statistical weights 
and consideration of other Lyman series transitions (Ly/ and above) balances revision in luminos- 
ity. Additional uncertainty surrounds the possibility of UV-absorbing dust in the early IGM due to 
early stellar processing and expolusion. 

A straight-forward calculation shows that existing instruments such as the GBT, Parkes, and 
the most compact configuration of the ATNF, could potentially detect, or set interesting limits on, 
the neutral fraction of the IGM at z ~ 6 in reasonable integration times ~ 100 hours. The primary 
challenge in detecting lsij2 — Ip^/i absorption will be achieving the very high spectral dynamic 
range required, of order 10^. 

4. Radio studies of the first galaxies 

4.1 Introduction 

The last few years has seen remarkable advances in the discovery of luminous objects into 
cosmic reionization, ie. z > 6. This includes over 100 star forming galaxies discovered through the 
Ly-break and Lya, and related, techniques (Ellis 2006), and some tens of QSOs (Fan et al. 2006). 
An important point to keep in mind when observing the first galaxies is that the GP effect precludes 
detection of these sources at observing wavelengths < 900nm. Hence, study of the first luminous 
objects is the regime of near-lR through radio astronopy, and hard X-rays. 

Radio astronomy has also detected objects at z 6 at cm and mm wavelengths in line and 
continuum emission. In this section 1 want to emphasize complementarity. Proper study of the 
complex process of galaxy formation requires observations across the electromagnetic spectrum. 
Near-IR observations reveal the ionized gas, stars, and AGN. Far IR through (sub)mm observations 
reveal the dust, higher order molecular Une emission, and low excitation fine structure line emis- 
sion. The higher cm frequencies reveal the lower order molecular transitions, which provide the 
best estimates of total gas mass. And at longer cm wavelengths one observes synchrotron emis- 
sion, which is a dust-unbiased measure of the star formation rate and distribution, and of obscured 
AGN. The mm and cm studies are crucial to study the cool thermal material, the fuel for galaxy 
formation, and to probe galaxies in their earliest, dust-enshrouded stages of formation. 

In this section I also emphasize the fact that a full SKA is not absolutely required for this 
program, and that a 10% Phase I SKA operating up to ~ 40 GHz would be a major step toward 
meeting this KSP (Carilli 2006). 

4.2 The case of 1148+5251 at z = 6.42 

To illustrate the power of our radio studies, we describe in some detail the particularly en- 
lightening example of the most distant SDSS quasar, Jl 148+5251 at z = 6.419, a luminous broad 
Une AGN with a black hole mass ~ 10^ M© (Fan et al. 2004). The host galaxy has been detected 
in thermal dust, and non-thermal radio continuum, as well as molecular line emission, including 



15 



Cosmic reionization and first light 



C.L. Carilli 




Cn ISSim IRAMPdBI 



Figure 10: Images of the most distant SDSS QSO Jl 148+5251 at z=6.42. Left: The optical/nearlR SDSS 
color image (Fan et al. 2004). Center: The CO 3-2 emission from the host galaxy observed with NRAO's 
Very Large Array (Walter et al. 2004). Right: The Cll 158um emission observed with IRAM's Plateau 
de Bure Interferometer (Walter et al. in prep.). This supermassive black hole (black hole mass le9 solar 
masses) is enveloped in a giant molecular gas cloud, with a molecular gas mass 2el0 solar masses, dis- 
tributed over a physical scale of about 5.5kpc. The extended Cll emission suggests active star formation 
across this region, with a total star formation rate 3000 solar masses per year (Maiolino et al. 2006). 



multiple CO transitions, and interesting limits to dense molecular gas tracers, such as HCN (Walter 
et al. 2003; Bertoldi et al. 2003; Carilli et al. 2004; Riechers et al. 2007). 

Most recently, we have detected [CII] fine structure line emission from Jl 148-1-5251 (Maiolino 
et al. 2005). Fine structure lines are the dominant cooling mechanism for interstellar gas, and hence 
a key diagnostic of the energetics of the interstellar medium. Figure 10 shows the high resolution 
image of the [CII] and CO emission from J1148-I-5251. The [CII] and CO emission are co-spatial, 
and extended over a scale ~ 5.5 kpc, suggesting distributed gas heating, and hence star formation, 
on kpc-scales. 

Figure 1 1 shows the rest-frame radio through near-IR SED of Jl 148-1-5251 (Wang et al. 2007). 
The source shows a clear excess in the FIR over the expected SED of a typical low z, optically 
selected QSO. This SED from the FIR through the radio is consistent with that expected for an 
active star forming galaxy with a dust temperature of ~ 50K, and follows the radio-FIR correlation 
for star forming galaxies. Figure 4 also shows the CO excitation. The gas excitation is reasonable 
fit by an LVG model implying dense (> 10^ cm^-'), warm (> 50K) molecular gas, comparable to 
what is seen in the centers of nearby nuclear starburst galaxies (Bertoldi et al. 2004). 

These observations of Jl 148-1-5251, and other z ~ 6 QSOs (CariUi et al. 2007; Wang et al. 
2007), demonstrate that large reservoirs of dust and metal enriched atomic and molecular gas can 
exist in massive galaxies within 1 Gyr of the Big Bang. The current observations suggest active 
star formation in the quasar host galaxy, with a massive star formation rate of order 10^ M0 year^^ 
adequate to form a large elliptical galaxy in a dynamical timescale 10^ years. 

These observations of the host galaxies of the most distant quasars have inspired a number of 
major theoretical efforts aimed at delineating the formation of massive galaxies and supermassive 
blacks at the earliest epochs. In particular, the theory groups at Harvard and Arizona have per- 
formed state-of-the-art numerical simulations which follow the formation of the largest galaxies 
in the early Universe. In a series of recent papers (eg. Li et al. 2007; Narayanan et al. 2007), 
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Figure 11: Left: The radio through near-IR SED of J 1148+5251 (Wang et al. 2007). Right: The CO 
excitation in Jl 148+5215 (Bertoldi et al. 2003). 

they show that early massive galaxy and SMBH formation is possible in rare peaks in the cosmic 
density field, through a series of gas-rich, massive mergers starting at z 14. These systems evolve 
into massive galaxies at the centers of the most massive cluster environments seen today (~ 10^^ 
M0). These results are generally consistent with the idea of 'downsizing' in both galaxy and su- 
permassive black hole formation (Cowie et al. 1996; Heckman et al. 2004), meaning that the most 
massive black holes (> 10^ M©) and galaxies (> 10^^ M0) may form at high redshift in extreme, 
gas rich mergers/accretion events, in rare peaks in the cosmic density field. 

4.3 The Coming Revolution: ALMA and a Phase I SKA 

While radio astronomy has pushed back into cosmic reionization through these molecular line 
and continuum observations, current studies are at the limit of the capabilities of existing millimeter 
and radio telescopes. Even at these limits, these experiments stretch current instrumentation to the 
extreme limit, such that only rare and pathologic objects are detectable, ie. ultra- to hyper-luminous 
IR galaxies, with star formation rates > 100 M© year^^ Fortunately, in the coming few years there 
will be dramatic improvements in observational capabilities at submillimeter through centimeter 
wavelengths that should revolutionize studies of the earliest galaxies. 

An order of magnitude, or more, increase in collecting area at mm and cm wavelengths is 
required to enable the study of the first 'normal' galaxies within the epoch of reionization (EoR), 
eg. the z > 6 Lya galaxies currently being discovered in deep near-IR surveys (SFRs ~ 10 to 100 
M0 year^'; eg. Murayama et al. 2007). At mm wavelengths this increase in capability will be 
realized when the Atacama Large Millimeter Array becomes operational early in the next decade. 
ALMA represents close to two orders of magnitude improvement in all areas of submillimeter 
astronomy, including: spatial resolution, spectral capabilities, and sensitivity. At cm wavelengths, 
the Expanded Very Large Array opens up the full frequency range from 1 to 50 GHz with order of 
magnitude improved continuum sensitivity, and dramatically improved spectral search capabilities. 
Unfortunately, the EVLA will still have the same total collecting area as the VLA. An increase in 
collecting area by roughly an order of magnitude over the EVLA is required at cm wavelengths to 
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match advances at other wavelength ranges, and perform the key, complementary science afford at 
the short cm wavelengths. 

The key areas of science explored by these new facihties includes: 

• Dust and molecular gas: ALMA, and a Phase I SKA (10% demonstrator ~ lOx EVLA 
collecting area) will reveal the dust and molecular gas in samples of normal star forming 
galaxies at z > 6, including Ly alpha emitters and Ly -break galaxies. These observations 
constrain the timescale for metal and dust enrichment in the most distant galaxies, and probe 
a key element in galaxy formation, namely the fundamental fuel for star formation. 

• Interstellar medium physics: The physical conditions in the ISM of these galaxies (density, 
temperature, abundances) can be studied in detail, including gas excitation conditions, dust 
spectral energy distributions, and observations of dense gas tracers such as HCN and HCO+, 
directly associated with star forming regions. 

• Fine structure lines: The [Cll] 158/xm line, and other atomic fine structure lines, are the 
dominant ISM gas cooling fine structure lines. These lines have great promise for exploring 
star formation, ISM physics, and galaxy dynamics, in the first normal galaxies. I expect these 
Unes to be the 'work-horse' lines for studies of the first normal galaxies using ALMA. 

• Galaxy dynamics: High resolution imaging of the distribution of molecular and ionized gas 
is the most direct method to determine the gas dynamics and dark matter content of the 
earliest galaxies. These observations also present the only method with which to test the 
M-<7 relation for galaxies and black holes out to the highest redshifts. 

• Star formation: The centimeter continuum emission associated with star formation in these 
galaxies can be observed, including millarcsecond resolution imaging using Very Long Base- 
line Interferometry. 

4.4 Continuum examples 

I illustrate the continuum studies with two recent examples. First, Figure 12 shows the results 
of a median stacking analysis of a sample of 6500 z ~ 3 LBGs in the COSMOS field. We obtain 
a robust statistical detection at 1.4 GHz of 0.9 ±0.15 /iJy. Through the stacking analysis we are 
reaching sub-/xJy sensitivity, ie. SKA-like depths, on galaxy populations at high redshift. The 
interesting implication is that the median star formation rate derived from the radio equals that 
derived from the median rest-frame UV emission, assuming the standard factor five dust correction 
for the UV emission. In other words, these data provide a critical, independent check on the 
standard reddening factor for high z LBGs (Carilh et al. in prep). We have also identified the first 
z > 4 submm galaxy, or massive dusty starburst, in these LBG samples (Capak et al. in prep). 

Second, Figure 12 also shows a VLBA image with 28 /iJy sensitivity of the most distant radio 
loud AGN, the z = 6.12 QSO J1427-I-3312 (McGreer et al. 2008; Momjian et al. 2008, in prep). The 
source has a steep spectral index of 4 = — 1.1. The source is comprised of two dominant com- 
ponents separated by ~ 170pc, with a flux density ratio of about 4/1. Both components are clearly 
resolved by these observations, with sizes ~ 60pc. These data are consistent with J1427-I-3312 be- 
ing a Compact Symmetric Object. CSOs are thought to be very young radio sources, < 10^ years. 
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Figure 12: Left: A median stacked image at 1.4 GHz of 6500 Lyman Break galaxies from the COSMOS 
field (z 3, U drop-outs; Carilli et al. in prep). The contour levels are: -0.75, -0.5, -0.25, 0.25, 0.5, 0.75 jxiy 
beam"i. Right: A VLBA image at 1.4 GHz resolution of the most distant radio loud AGN, J1427-H3312 at 
z = 6.12. The contour levels are: -3,3,6,9,12 times the rms noise of 28/1 Jy beam~^ (Momjian et al. in prep). 



highly confined by the dense ISM of the host galaxy (Taylor et al 1996; Readhead et al. 1996). The 
conversion efficiency of jet kinetic energy to radio luminosity is thought to increase in such highly 
confined sources. 

While it is dangerous to draw general conclusions from a single source, we can speculate that 
at the highest redshifts a higher fraction of radio jets will be CSOs, due to the higher gas densities 
in the earliest galaxies. Kinetic feedback from such jets can affect significantly the evolution of 
the host galaxy, perhaps inducing star formation, and eventually driving enriched material into the 
IGM. These radio source also present an ideal opportunity for observations of HI 21cm absorption 
during cosmic reionization (Carilli et al. 2002). 

I should also point out that VLBA observations at 1 .4 GHz are reaching (rest frame) brightness 
temperature sensitivities at mas resolution approaching 10^ K (Momjian et al. 2007). This limit 
is critical, since it is the dividing point between the brightness temperature expected for an AGN 
versus that predicted for a starburst (Condon et al. 1991). 

4.5 Complementarity and sensitivity 

There are a few important factors to keep in mind when considering SKA phase I studies of 
high z molecular line emission, in comparison to other, nearer term, telescopes: 

• When considering the very first galaxies, into the EoR (z > 6), mm telescopes such as ALMA 
will be limited to studying very high order transitions, eg. at z = 6 the lowest ALMA band 
(80 - 100 GHz) redshifts to 690 GHz, corresponding to CO 6-5, or HCN 7-6. It is unclear 
that these higher order transitions will be excited in an average galaxy, in particular for 
the high dipole moment molecules such as HCN, since the required densities for excitation 
become extreme (>> 10^ cm~^).-^ For comparison, in normal galaxies, such as the Milky 

■^The CMB will also contribute to excitation. For example, at z=6, T(CMB) = 19K, while the typical CO excitation 
temperature in an active star forming galaxy is ~ 30 to 40K. 
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Way, constant Tg holds for the integrated CO emission roughly up to CO 3-2, above which 
the lines are sub-thermally excited. 

• The SKA phase I design corresponds to a collecting area roughly an order of magnitude 
larger than the EVLA. While the EVLA continuum sensitivity at higher frequencies is a fac- 
tor 9 improvement over the existing VLA due to the increased bandwidth, the EVLA will 
have the same collecting area as the VLA, and hence will have at best marginally improved 
sensitivity for spectral line (ie. band-limited) studies, perhaps 50% due to small improve- 
ments in receivers and antenna optics. What the EVLA does provide for line science is 
improved search capabilities, in terms of new receiver bands and much wider bandwidths 
with many more channels. 

• It should be emphasized that the natural high frequency limit for cm telescopes is ~ 45 GHz, 
as set by the atmospheric O2 line. If the SKA phase I was extended to this natural hmit, 
the effective sensitivity is increased by another factor four for thermal objects (assuming 
constant brightness temperature) relative to a 22 GHz design, improving the SKA phase I 
to a 40% SKA. A higher frequency would also decrease the necessary longest basehne by a 
factor two at fixed resolution, thereby reducing the costs of array configuration, connectivity, 
and data transmission. 

• Beyond pointed observations of known high-z objects, the much wider field of view and 
fractional bandwidth of the SKA relative to eg. ALMA, will allow for efficient surveys for 
molecular emission line galaxies from z ~ 2 into the reionization (Carilli & Blain 2003). 

The key point is that, for the first galaxies, the fine structure lines, such as [CII] 158;Um will 
be the work-horse lines for ALMA, while studies of molecular gas, the fundamental fuel for star 
formation in galaxies, will be the regime of short cm wavelength telescopes. 

Figure 13 shows the sensitivity to molecular line emission for the EVLA, the SKA phase I, and 
ALMA. Also shown is a model for the expected continuum and CO line emission from an active 
star forming galaxy with a total IR luminosity ~ 10% that of Arp 220 at z = 5, corresponding to 
a star formation rate of a few lO's M© year~^ Such a galaxy would be comparable to the Lya 
galaxies seen in deep optical surveys at z 6, ie. characteristic of the high z galaxy population 
(Murayama et al. 2007). I have assumed a CO excitation ladder similar to that seen in the few high 
z molecular line emitting galaxies currently known (this may be biased to higher excitation). 

Figure 1 shows that the SKA phase I will be able to detect the CO 1-0 and 2-1 emission from 
'normal' star forming galaxies out to very high redshifts, into the epoch of cosmic reionization. 
ALMA will easily detect the higher order transitions, in sources where they are excited. ALMA 
will also be able to study the low excitation fine structure line emission from primeval galaxies 
([CII] 158Aim, [01] 63Aim). 

For completeness, in Figure 14, I show the continuum sensitivities of these telescopes. The 
high frequency SKA phase I will easily detect the continuum emission over the rest frame fre- 
quency range of 50 to 200 GHz range. The lower end of this range is where free-free emission will 
dominate, providing the most direct, dust-unbiased, measure of the star formation rate. And the 
upper end corresponds to the cold dust emission. 
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Figure 13: The dotted curves show the (Id) sensitivity to spectral hne emission for the VLA, SKA phase I, 
and ALMA in 12 hours. The solid curve is the emission spectrum, including continuum and CO lines, for a 
source with the luminosity of Arp 220 (ie. a source with FIR - 10'2LQ)atz = 5. 
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Figure 14: The solid curves show the (la) sensitivity to continuum emission for the EVLA, SKA phase I, 
ALMA, Spitzer, and the JWST, in 12 hours. The dashed curves are the continuum emission spectrum for a 
source with the luminosity of Ai-p 220 (ie. a source with FIR ~ lO'^ L©) atz — 2, 5, and 8. 
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Overall, the EVLA Phase ru and ALMA provide the necessary complement to large ground 
and space-based optical and IR telescopes, enabling a panchromatic study of the formation of 
normal galaxies back to the EoR. 
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